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ABSTRACT 

In this work we estimate the accretion-disk mass in the specific scenario of binary-neutron-star- 
merger with current observational data. Assuming that the outflows of short Gamma-ray Bursts 
(GRBs) are driven via neutrino-antineutrino annihilation we estimate the disk mass of about half of 
short bursts in the sample to be ~ 0.01 — 0.1 Mq, in agreement with that obtained in the numerical 
simulations. Massive disks (~ several 0.1 Mq) found in some other short GRBs may point to the more 
efflcient magnetic process of energy extraction or the neutron star and black hole binary progenitor. 
Our results suggest that some short bursts may be really due to the coalescence of double neutron stars 
and are promising gravitational wave radiation sources. For future short GRBs with simultaneous 
gravitational-wave detections, the disk mass may be reliably inferred and the validity of our approach 
will be tested. We also propose a method to constrain the initial radius of a baryonic outflow where 
it is launched (Ro) without the need of identifying an ideal thermal spectrum. We then apply it to 
GRB 090510 and get that Rq < 6.5 x 10^(rph/2000)~'* cm, suggesting that the central engine is a 
black hole with a mass < 22 MQ(rph/2000)~^, where Fph is the bulk Lorentz factor of the outflow at 
the photospheric radius. 

Subject headings: Gamma rays: general — Radiation mechanisms: non-thermal — Accretion, accretion 
disk 



1. INTRODUCTION 

Since the discovery of the afterglow and then the 
measurement of the redshift in 1997, our understanding 
of Gamma-ray Bursts (GRBs) has been r evolutionized 
(see lPiran1l2004l : [Zhang fc Meszaros 112004 for reviews). 
As usual some aspects are understood better than others. 
For example the late time {t > 10^ s) afterglow emission 
is likely dominated by the radiation of the electrons ac- 
celerated by external forward shock while the earlier af- 
terglow emission may consist of at least two components, 
including that powered by the prolonged activity of the 
central engine and the external shock emission (e.g. , 
Fan & Wei |j2005l : iMeszaros] [20061 : IZhang et all [20061 : 
Zhang II2007| ). The origin of the prompt 7— ray emis- 



sion is less clear. Widely discussed scenarios include the 
internal shock model, the internal magnetic energy dissi- 
pation models, and the photospheric models. The GRBs' 
central engine is either a stellar black hole surrounded 
by a hyper-accreting disk (iMacFadvcn fc Woo slev l[l999| ) 
or a quickly rotating magnetar (Usov 1992). For a 
magnetar-like central engine the outflow is expected to 
be Poynting-flux dominated. In the case of a stellar 
black hole surrounded by a hyper-accreting disk, the 
outflow could be either baryonic or Poynting-flux domi- 
nated, depending on the energy extr action process (see 
lPiranl[200l IZhang fc Meszaros1l2004 for reviews). Such 
a kind of central engine is usually characterized by the 
mass and the spin of the black hole, and the rate of ac- 
cretion onto the black hole. Our current knowledge of 
these physical parameters is mainly from the numerical 



simulations since current electromagnetic data alone can 
not break the degeneracies between the parameters (see 
section 12. 1[) . For some specific long bursts (for which 
the duration is longer than ^ 2 s) the situation is bet- 
ter. The modeling of their associated supernovae sheds 
some light on the mass of the central remnant (e.g., 
iDeng et al. 1 120051) . However, for short GRBs (with a 
duration less than 2 s) no bright associated-supernova 
has been detected. In this work we investigate whether 
it is possible to make some progress with some spe- 
cific assumptions. Wc concentrate on the binary neu- 
tron star merger model which has been supported by the 
host galaxy observations and by the non-detection of ac- 
companying bright supe rnova for some short bursts (e.g., 
Gehrels'et al. 2005; Fo x et al. 1 120051: iBarthelmv et al. I 
2005;iBerger et al. l2005HHiorth et al. II2005D . The other 



possibility that so me short even ts might have a mas- 
sive star origin (e.g. .IZhang et al. 2003; Fan ct al. 2005|; 
IZhang et aLll2009l: IVirgih et al. 112011; ,Panaitescu„201ir) 
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won't be addressed. 

Besides the physical parameters of the central en- 
gine, the initial radius of the outflow where it is launched 
(Ro) is an important parameter revealing the physical 
process taking place at the center of the burster. In 
the coUapsar scenario, the interaction of the acceler- 
ated/cooled ejecta with the envelo pe material may give 
rise to a re-born hot fireball (e.g.. ILazzati et al. II2009I) 
and then the derived Rq marks the site of the interac- 
tion, which is much larger than the initial size of the 
ejecta. If such an interaction is ignorable, Rq imposes an 
independent though rough constraint on the mass of the 
central black hole. For a baryonic outflow the accelera- 
tion and the energy dissipation processes are well under- 
stood (Phan ct al. 1993; Meszaros ct al. 1993). With 
a reliable thermal component identified in the prompt 



spectrum, it is possible to c onstrain the bulk Lorentz 
factor and Rq of the shells (jPe'er et al. II2007D . Such 
a goal was achieved for some long GRBs, in particular 
GRB 090902 B. and a typical -Rn ~ 10^ - 10^ cm was 
inferred fe.g.. IPe'er et al. 1 120071: iRvde et al. 1 12010'). For 
short GRBs, no reliable thermal spectrum component 
has been identified. A unique candidate is GRB 090510 
characterized by a very soft Me V spectrum and a rather 
hard GeV radiation component ()Gao et al. Il2009f) . But 
the peculiar spectrum of this short burst, unlike the long 
event GRB 090902B, can not be reasonably fitted by a 
thermal component superposed by a power-law compo- 
nent (B. B. Zhang & B. Zhang, 2011 private communi- 
cation). Therefore one purpose of this work is to find a 
way to estimate Rq without the need of identifying an 
ideal thermal component. 

This work is structured as the following. In section 
2 we discuss the difficulty of estimating the physical pa- 
rameters of the central engine with current observational 
data, and then show that some progress is achievable in 
the specific double neutron star merger scenario. The 
mass of the accretion disk of some short GRBs has been 
estimated. In section 3 we present a method, in which 
an ideal thermal signature is not needed, to estimate Rq 
with the assumption that the outflow is baryonic. We 
apply the method to GRB 090510 and then constrain 
Rq. Our results are summarized in section 4 with some 
discussions. 

2. THE MASS OF THE ACCRETION DISK OF SOME 
SHORT GRBS 

2.1. The difficulty of estimating physical parameters of 
the central engine with current limited data 

One nascent stellar black hole surrounded by a 
hyper-accreting disk, the widely adopted central engine 
of GRBs, is characterized by some important parame- 
ters, including the mass and the spin of the central black 
hole (Mgjj and a, where a = J/GM'^^c, J is the angu- 
lar momentum of the black hole and c is the speed of 
light), the mass of accretion disk Afdisk or alternatively 
the accretion rate M. 

The accretion disk is so hot that the energy 
loss may be mainly through neutrino and anti- 
neutrino radiation. The later neutrino and anti- 
neutr ino annihilation m ay launch a baryonic fire- 
ball (jEichler et al. I [l989[ ). The annihilation luminos- 
ity has been extensively investigat ed in the litera- 
ture ( e.g.. iMacFadvenfc Wooslev l[l99 9: Popham et al. I 
[19991: iLiu et al. 1 120071: iZalamea fc BeloborodovilmTI T 
In general the luminosity is a function of M^^, a, M, 
and possibly also the vertical structure of the disk 0. In 
this work we adopt an empiric al relation proposed by 
IZalamea fc Beloborodov I ()201lD . which reads 

L.,«10^^ergs-ix-r(^)-^/^ 

^ The same holds for the models in which it is the mag- 
netic process rather than the neutrino process to extract 
the energy. For example in the Bland ford-Znajek process 
HBlandford &i Znaiek I [1977] ; [Lee et al. 1 12000| '|. the luminosity of 
the electromagnetic outflow can be estimated by L^^ ~ 2.5 X 

lO'*^ erg (a/0.5)2(MBjj/2.7Mgjj)^BH,i5' ™^*^'''^ ^a ~ ^'^ ^ 
lO^'^ Gauss (M/0.01 M© s-^)^^'^R~\ is the magnetic field 

strength on the horizon and R^ = (1 + Vl — a''^)r^/2. 



■ 0, for M < Mig„; 

771^/4, for Mign <M < Aftrap; 



(1) 



^ trap 7 



where the 



Mtrap 

The coefficients K\ 



accretion rate m = M/Mq s ^, 
r„,,(a)/rg, Mig„ = i^ig„(a/0.1)5/3, 
-f^trap(a/0.1)^/"^, and a is the viscosity, 
ign and Xtrap are functions of the 



black hole spin a. For a = 0, Kig^ = 0.071 Mq s ^ and 



K, 



trap 



9.3 M, 



0S 



For g = 0.95, Kjgn = 0.021 Mq s' 



and Ktre^p = 1.8 Mq s-1 ()Chen fc Beloborodov I [20071) 
The radius of last stable orbit r^^s is (jBardeen et al. I 

Tms = rg{3 + Z2T [(3 - Zi)(3 + Zi + 2Z2)]i/2|/2^ (2) 



where 



Zi = 1 + (1 - a2)i/3[(i + a)i/3 + (1 - a)'/% 



and 



Z2 = (3a^ + Zf) 



2U/2 



= 3rg, while for a = 1 we have 



For a = we have ' 

'"ms = rg/2 or r,ns = 9rg/2 (retrograde). The retrograde 
case is irrelevant to the case of accretion disk and will 
not be discussed in this work any longer. 

With the observational data the energy output of 
the central engine i?out and hence Li,p can be reasonably 
inferred, which however is not enough to break the degen- 
eracies among the parameters a, M^^ and to. As shown 
in footnote [31 the same applies to the magnetic process 
of energy extraction. That's why it is rather hard to con- 
strain the physical parameters of the central engine with 
current electromagnetic data and our knowledge of the 
central engine is mainly from the numerical simulations. 
Fortunately, the specific scenario of binary-neutron- star 
merger could be an exception, for which a preliminary 
probe is plausible. 

2.2. An exception: the double neutron star merger 
scenario 

In the specific binary-neutron-s tar merger model , 
also the lead ing one, for short GRBs (.Eichler et al. 1119891 : 
[Nakar 2007|; iLee fc Ramirez- Ruiz Il2007f ). the mass of the 
formed black hole and its spin parameter can be rela- 
tively reasonably evaluated. 

The mass of the central engine is expected to be 
close to the total mass of the progenitors since the 
mass ejection during t he merger is expected to be tiny 
(jRosswog et al. I ri999f ) and the mass range of Neutron 
Stars is relatively narrow. For the 10 neutron-star- 
binaries well studied so far, the total-mass of the b inaries 
ranges from 2.57Mq to 2.83M0 (Kiziltan et al. l[2010) . 
Therefore the mass of the nascent black hole formed in 
the merger is expected to be close to M^^ ^ 2.1 Mq. 

The spin parameter of the formed blac k hole can 
be semi-quantitatively estimated. Following iLee et aLl 
(2000) and for simplicity we assume that the double 
neutron stars have a similar mass Mj^g ^ 1.4 Mq, 
the orbital angular momentum of the binary system 
is ^binary = -^GAf ^^Tt /2, whcrc Tt is the tidal ra- 
dius. The spin parameter of the formed black hole 
can be expressed by a = a;yJbinary/[(2yAfNs)^G/c] = 



(a;/42/)^rt/rg 



where r„ 



2GAfpjs/c^, assuming 



4(l-J9)/(j9+2) 



that a fraction (xy) of the orbital angular momen- 
tum of the neutron star binary goes into the nascent 
black hole keeping a fraction y of the total mass. For 
the binary-neutron- star merger scenario, the tidal ra- 
dius is rt « 6rg_,^g ([Haensel et al. I [T991D and we have 
a = 0.61x/y, about 10% smaller than that suggested 
by iLee et al I (|2OOO0 . As already mentioned, y is close 
to 1. Let's estimate x. The gravitational radiation 
changes the energy of the binary system at a rate 
dE/dt « 1.9 X 10^^(rt/6rg^^g)^^ erg/s and the merger 
takes a time Ai ~ 1.4 ms (j't/Sr-g^^g)'* (jHaensel et al. I 
|1991|) . The corresponding change of the angular mo- 
mentum AJbinary 



{dE/dt)At/y/2GM^Jrl We then 



have X = 1 — AJbinary /-/binary '^ 0.9. Therefore we have 
a ^ 0.55 . i.e., the formed black hole rotates rapidly. As 
found in iLee et al. I ()2000D . for the double neutron stars 
having different masses M^g_-^ and Mj^g_2, roughly one 

has a oc 7^ ^ M^,_,M^/lJiM,,_, + M,,_J"/6. For 
the 10 neutron-star- binaries discussed in lKiziltan et al.1 
(|2Q10D . one finds out that TZ ranges from 0.275 to 0.282, 
i.e., a is insensitive to the mass ratio of the binaries. The 
above semi-quantitative analysis is in agreement with the 
numerical simulation, in which people found a typical 
a ~ 0.78, weakly depending on the total mass and the 
mass ratio of the binary neutron stars (jKiuchi et al. I 
[20091) . 



For a = 0.78, we have Vms 

(for Afign <M < Mtrap) 



1.45r, 



g' 



= 1.45 and 



-l51 



,9/4 



1.45 



r'H^r"'. (3) 



A typical Mdisk ~ 10"''-10"^ Mq is found in the nu- 
mer ical simulation of th e binary-neutron-star coalescence 
(see lKiuchi et al. 1120091 and the references therein). The 
binary-neutron-star merger hypothesis will be supported 
if our Afdisk estimated with the observational data is 
within such a mass range. In the foreseeable future the 
binary-neutron-star merger model for short GRBs could 
be directly tested by the Gravitational wave data and 
the mass of the binaries as well as the formed disk could 
be inferred (jKiuchi et al. |[2010l ). Therefore the validity 
of our following approach will be directly tested. 

2.3. A simple approach and case studies 

A simple approach. Below we take the simplest 
approach to estimate L^p. The isotropic-equivalent ki- 
netic energy of the outflow powering long-lasting after- 
glow (i?k,iso) and the opening angle of the ejecta 9j can 
be derived fro m the mo deling of the multi-wavelength 
afterglow data (jPanaite scu 2006). However, in a good 
fraction of short GRBs, such a goal is not achievable due 
to the lack of prompt observations. Fortunately, for the 
X-ray emission above both the typical synchrotron ra- 
diation frequency and the cooling frequency, the flux is 
independent of the poorly constrained number density of 
the medi um n and the X-ray lumin osity is a good probe 
of Sk.iso (jKumar 112001 . Following iFan fc Piran I ()2006[ ) 
we take 



-Ek 



lO^-" erg £ 



4/(p+2) ^ 1 + ^ ^ ,-(p-2)/(p+2) 
X,46 



(^^) 



(i + r)4/(P+2)^ 



(4) 



where Cx is the X-ray afterglow luminosity at t = 10 
hours after the trigger of the burst, e^ and ee are the 
fractions of shock energy given to the electrons and 
magnetic field respectively, Y is the Compton parame- 
ter, and p ~ 2 is the energy distribution index of the 
shock-accelerated electrons and is constrained by the X- 
ray spectrum. The convention Qn ~ Q/10" has been 
adopted here and throughout this work except for some 
specific notations 

The jet opening angle is estimated to be (jFrail et al. I 

[2OOII) 



0.076(tj/lday)3/^[(l-Hz)/2]- 



-3/8^-1/8 „l/8 



'E, 



k,iso,5l' 



(5) 



where tj is the jet break time. We then estimate the 
"intrinsic" power released by the GRB central engine as 



E^ 



(1 + z)iE^,iso + Sk,iso)ef /(2ract), 



(6) 



where Tact is the duration of the activity of the cen- 
tral engine. In reality, £^out is just a fraction {J-' < 0.3) 
of the total neutrino-antineutrino annihilation luminos- 
ity outside the horizon of the rotating black hole, i.e., 
£^out = J^Li,i, (e.g.. lAlov et al. II2005D . The real duration 
of the main activity of the central engine might be shorter 
than the duration of the prompt emission Tgo by a factor 
7?. > 1, i.e.. Tact = Tgo/TZ, because different propagation 
velocities of the front and rears ends will lead to a radial 
stretc hing of th e ultra- relativistic ejecta (see section 4.1 
of Alo v et al. I '2005. for more details). Hence we have 

(for Afign <M < A/trap) 



• ncodf r^(l + ^)(-^7,iso,51 +-Ek,iso,5l)6'? 4/9 

mwO.53 Mq [ ^1 ' 



^To, 



1.45'^ ^2.7Mq^ ' 



(7) 



and 



A'/disk«O.53Af0 [ 



(^7,iso,51 + £'k,iso,5l)^j 



|4/9 



J" 



^90 ]5/9|- ^ms x2.1j^ ^"BH ^^2/3 



2.1. Mbh .2/ 



'{l + z)n' '1.45' '2.7AfG 



(8) 



Since 7?. > 1 and J- < 1, their impacts on estimating 
Afdisk are partly canceled. For a — 0.6, we have 
Xms ~ 1-9, Afdisk given in eq.(j8]) will be enhanced by a 
factor of ^ 1.8. 

Case studies. So far we have 10 short GRBs, 
as listed in Tab.l, having relatively abundant afterglow 
data, with which we can estimate -Eout and then Afdisk- 
For simplicity, instead of discussing all bursts one by one, 
below we focus on a few special events. 

GRB 051 221 A, a burst with a duration 2 = 1.4 s, 
was at a redshift z — 0.5465 and is distinguished by a 
long-lasting X-ray flat segment in the afterglow. Such 
a flat segment could be due t o either the energy in- 
jection from the cent ral engine (jSoderberg et al. 1120061 : 
iBurrows et al. I [2006f) or the forward shock emission of 
the emerging wid e-component of the two component jet 
pin et al. II2007D . The afterglow parameters reported in 



the literature are different. In the analysis we take the 
parameters obtained in the two-component jet model- 
ing, in which the long activity of the central engine is 
not needed. If we take the so mewhat more conserva- 
tive es timate £^01 « 3 x 10^^ erg (Sodcrbe rg et al. II2006I : 
[Burrows ct al. 2006) , the mass of the accretion disk will 
be reduced by a factor of 0.6. 

GRB 090510, a burst at a redshfit z = 0.903, is the 
most energetic short event ever recorded and is also re- 
markable for its long-lasting GcV emission that is likely 
powered by the external f o rward shock (Gao et al. 2009; 
IDe Pasauale et al. I l2010t ICorsi et al. I I2OIOI ) . The self- 
consistent interpretation of the GeV/X-ray/optical data 
is not an easy task and the parameters are found to 
be somewhat unusual. M t '^ 1200 s after the trig- 
ger of the burst, the X-ray and the optical afterglow 
emission change the dec line behaviors achromatically 
(|De Pasquale et al~ll2010t ). Such changes are most likely 
due to the jet effect, that is, the edge of the outflow en- 
ters our line of sight and the visible emitting region can 
not be approximated as a spherical surface any longer. 
The jet opening angle is as small a s ~ .006 (jGao et al. I 
|2009; Corsi ct al. 2010; He et an i2011| ). which is about 
one order of magnitude smaller than that of other bursts 
see Table 1) or th at found in the numerical simulation 
Alov et al~ll2005[ ). It is unclear how such a narrow col- 
limation is reached. Nevertheless a similar narrow col- 
limation was identified in the a fterglow modeling of the 
naked-eye burst GRB 0803 19B (Racusin et al J 120081). 

GRB 1008 16 A, as a burst at z == 0.8035, the local 
duration is 2.8/(1 + z) = 1.55 s, consistent with being 
a short burst. The result of the Swift-BAT spectral 
lag analysis is also consistent with being a short hard 
burst, but the erro r bars are too large to be definitive 
(jOates et al. I l2010f ). With the X-ray light curve we 
take a jet break time ij > 2 x 10^ s. The most valuable 
information inferred from the afterglow data of this 
burst is likely the density profile of the circum-burst 
medium. The preliminary white band flux decline is 
~ t"ii5 for 100 s < i < 10^ s, steeper than the s i- 
multaneous t^^— like X-ray decline ()Oates et al. |[20Tol ). 
The spectral index of the X-ray afterglow photons is 
~ 1.03to;^g, suggesting that the X-ray emission is above 
the cooling frequency of the forward shock electrons 
for p ~ 2. If the slow-cooling fireball was expanding 
into the ISM-like medium and the typical synchrotron 
radiation frequency is below the observer's band, the 
optical afterglow emission should decline with the time 
as i""-^^— like, shallower than the X-ray decline, which 
is at odds with the data. For a free-wind medium, 
the optical dec line should be t~^'^ ^ , stee per than the 
X-ray decline ([Zhang fc Meszaros I |2004[ ). Therefore 
the current data favor the free-wind medium model, in 
which the progenitor should be a massive star rather 
than a pair of compact objects. If our speculation 
could be confirmed by the careful analysis of available 
afterglow data reported in GCNs, the collapsar origin 
of GRB 100816A with an intrinsic duration ^ 1.4 s 
would be established. In turn such a result would be 
in support of the hypothesis tha t collapsar can produce 
short events (jZhang et al. II2003D and the progenitors of 
short GRBs are diverse (iFan et al. 1 120051 : iZhang et al.l 



of bursts is likely invalid. 

As shown in Tab.l, for about half of short bursts 
in the sample, the accretion disk has a mass ^ 0.01 — 
0.1 Mq, well consistent with that found in the numer- 
ical simulations of double neutron star merger (e.g., 
iRosswogetaTI 120031: iKiuchi et al"! I2OO90 . For some 
other events, such as GRB 051221A and GRB 050724, 
the inferred Mdisk ~ several 0.1 Mq may be a bit mas- 
sive to form. This puzzle can be solved in either of the 
following scenarios. One is that the outflows of these 
short GRBs were launched via some more efficient mag- 
netic processes rather than the n eutrino mechanism, as 
speculated in the literature ie.s... iRosswog et al. 1 120031 : 
iLee fc Ramirez-Ruiz I [2OO7I ). With footnote H for a k. 
0.78, M^H = 2.7M0 and M = (0.1, 1.0)1/© s^^ we have 
L32 ^ (100, 10)Li.p. Consequently the disks with Afdisk 
about 10 or more times smaller than that presented in 
Tab.l may be enough to power these short events. The 
outflow launched in this way is Poynting-flux-dominated 
and the prompt emission due to the magnetic energy 
dissipation should have a high linear polarization degree. 
The other is that some short GRBs are from the neutron 
star-black hole merger for which a massive disk is possi- 
ble. Another possibility that can not be ruled out is that 
some short events might have a massive star orign. Since 
our estimated Afdisk is close to that found in the numer- 
ical simulations, the compact object merger scenario is 
likely viable, implying that some (possibly a consider- 
able fraction of) short GRBs may be really driven by the 
coalescence of double neutron stars and are promising 
gravitational wave radiation sources. 

3. ESTIMATING THE INITIAL RADIUS OF THE OUTFLOW 
3.1. The method 

The fol l owing approach is partly motivated by 
iPe'er et al. I ()2007[ ). In their work the photospheric ra- 
diation taking place in the radiation-dominated phase 
and matter-dominated phase have been investigated sep- 
arately while we solve the problems jointly. Furthermore 
we focus on constraining Rq in the absence of an ideal 
thermal signature, different from what did in the litera- 
ture. 

For a baryonic outflow, most of the initial ther- 
mal energy may have been converted into the kinetic 
energy of the baryons at the end of the acceleration 
(jShemi fc Piran 1 119901 ) but a (quasi-)ther nial em ission 
component is likely inevitable (|Paczvhsknil990() . The 
(quasi-)thermal emission is mainly from the photosphere 
at a radius i?ph which satisfies the following relations. 

B ased on its defini t ion, i?nh can be expres s ed as 
(e.g., j Paczviiskll [19901 : iDaigne fc Mochkovitch I [200I 
iJinet ah , 20ldr 



i?ph«4.5 X 10" cm L^iV-l^^^ 



i3.7x 10^X54/"^"^^ 



(9) 



l2009f ) . If so the calculation made in Table 1 on this kind 



where / = 3rph/477, rj is the initial dimensionless entropy, 
and Fph is the bulk Lorentz factor of the outfiow at the 
photospheric radius. 

T he photospher ic radius is also related to R n. Fo l- 
lowing [Piran et al.~[ ([1993D and [Meszaros et al. [ (|1993D , 
we introduce a, Ro, at which the bulk Lorentz factor of 



the outflow is To '-^ a few. With the parameter 

the acceleration calculation yields (|Piran et al. II1993I) 

Rph = RoiTo/Tpiy/^D^/^, (11) 

where r]o — eJ,/nJ,TOpC^ « v/^o- We then have -p = 
To/Fph + 3r2/4ryrph - 3r„/477 = 3(1 - /)r„/4/r, + 
9r^/16T]^f, the first term will be dominant as long as 
1 — / > 3ro/4?7, which is usually satisfied. So we have 

l/D^3il-f)To/Afv, OT<,/77«4//[3(l-/)], (12) 

with which we get 



1.3x10'" cm 7^3 i?o,7- 



(l_/)3/2' 
(13) 

where the relation Rq « Ro/^o has been used. 

Finally the photospheric radius is constrained by the 
observational data. Suppose the (quasi-)thermal emis- 
sion has a temperature Tobs and a fiux Fbb, we have 
47rr2j^ i?2^crTj^i^'' = ibb = 47ri:)^Fbb, which can be sim- 
plified as 

i?ph « [Fbb/aT^J 1/2(1 + z)-2rphi?L 



:1.3x 10^' cm K 



n _ .1/2 r(l±f)Zkil-2_^^3^^^^4) 



bb.-4L 



1 MeV 



where Tobs ~ ^p\iT'y-^/{1 + z) has been taken into account 
and and a is the Stefan-Boltzmann constant. 

Denoting the co-moving thermal energy density and 
the number density of the outflow at i?ph as e' j^ and 
n' -^ respectively, we have the (quasi-)thermal luminosity 



^bb 



<=ph 



4eph/3+"U™p'; 

n'^/D^ (jPiran et al. Ill993[ ). we have 



L. Since e'j^ — e'^jD^ and n' ^^ 



ibb 



e^ 



-,L 



L 



Ae'J'd + Dn'^TTipC^ 4/3 + DToh 



3(W)^ 
4 
(15) 



With Ibb = L\,\y/L < 3/4, the above relations yield / 
l-4ybb/3 and 77 = 3rph/[4(l-4Ybb/3)]. Combing eq.®, 
eq. p^ and cq. p^ we have 



Ro « 1.5 X lO'^ cm FbY,-4^bb 



1/2 ^3/2r(l+fy?obs, ; 



1 MeV 



-Cl,28: 

(16) 



ph 



lO'^ 



3 fy-l- 4/3^1/4 pl/8 r (l + g)Tobs .i/2pl/4 



1 MeV 



'L.28- 

(17) 

If the baryon loading is so low that at Rph the outflow 
is still radiation-dominated, most of the initial energy of 
the outflow will be lost via the thermal radiation and 
^bb '~ 3/4, for which Fph can not be reliably inferred, 
reflecting the well-established fact that both the observed 
temperature and the thermal radiation luminosity are 
constant until most of the initial energy of the outflow has 
been transferred into the kinetic energy o f the particles 
(jPiran et al. Ill993t iMeszaros et al. II1993D . For Fbb < 1 
(i.e., i?ph is far above the coasting radius ~ t]Ro), ea.pT)) 



sugges ts Fph oc V(jjj , in agreement with iPe'er et al. I 
((20071) . 

The above two equations flnally give 

Ro « 1.5 X IQS cm F^^f 3(^1^,1 - imY^l'Ph^^^.Dl^,. 

(18) 
Therefore, if Fph ami Ybb are obtainable with the observa- 
tional data, eg. il8\) provides us an independent estimate 
of Ro without the need of identifying a thermal compo- 
nent in the prompt spectrum. This is helpful since the 
physical processes taking place at i? < i?ph may be able 
to shape the thermal spectrum so significantly that the 
identificati on of an ideal therm al component is very dif- 
ficult (e.g.. iBeloborodov lIMol ). In fact, for short GRBs, 
no reliable thermal component has been identified so far. 
The disadvantage of our approach is that one can only 
get the time-averaged constraint. 

3.2. Application to GRB 090510 

GRB 090510 was detected by Fermi 7— ray tele- 
scope and Swift satellite simultaneously. Though the 
physical origin of the prompt GeV emission is not 
clear yet, the long lasting GeV afterglow emission 
is most likely the synchrotron rad i ation of the ex- 
ternal shock ( e.g., JGao et al. 1 12 009": 'G hirlanda et al. I 



20101: iCorsi et ah I I2010t iDeP asguale et all 120101: 



Kumar &: Barniol Duran II2010D . Supposing the prompt 



GeV and MeV emission are from the same region, the 
observation of GeV photons sets an upper limit on the 
optical depth for pair production and then suggests 
a bulk Lorentz factor of the emitting region > 1200 
(|De Pasauale et al. I [20lot lAbdo et aU [200I . To in- 
terpret the GeV emission at i > 2 s as the forward 
shock emission, a higher initial bulk Lorentz factor of 
the outfiow Pint > 1900 is required (IHe et al. 1 120111: 
IGhirlanda et~an mM ICorsi et al 1 12010|). Obviously 
Fph should be larger than Pint (it is straightforward to 
show that the outfiow shells with higher bulk Lorentz 
factor contribute more to the thermal emission). In the 
following estimate we take Fph '~ Pint ~ 2000. The time- 
averaged spectrum of GRB 090510 in the time interval 
0.5 — 1.0 s can be nicely fitted by a Band function plus 
a power-law component. The Band function component 
has an isotropic-equivalent energy £'Band,iso ~ 7 x 10^^ 
erg while the very hard power-law component (with a 
spectrum F^, oc i/~'^-^'^ in the energy range 10 keV— 10 
GeV) has an isotropic-equivalent energy i<^pL,iso ~ 5 x 
10^2 erg ()Abdo et al ] |2009 |). The afterglow modeling 
gives ^k,iso > 5 X 10^3 gj.g (iGao et al. I(200l IHe et al. I 
120111) . Clearly, only the Band function component may 
be relevant to the quasi-thermal radiation of the outflow. 
Hence we have a (quasi-)thermal radiation efficiency 

^bb < -£'Bandaso/(^kaso + ^PL,iso + -^'Bandaso) ~ 0.1 and 

the thermal radiation fiux Lbb ^ 6 x 10~^ erg s"^ cm^^. 
Substituting these values into ea.(|18p. we have 



Ro < 6.5 X 10^ cm r^l,,,Y^^ 



1/2 



(19) 



which is about 2 or more orders of magnitude smaller 
than that reported in the literature (jPe'er et al. 1 120071 : 
iRvde et al. I |2010[ !FI. suggesting that most energy has 

"* For one particular burst GRB 090902B, people found a promi- 
nent thermal signature in the prompt spectrum. The modeling 



6 



been deposited in a small cavity surrounding the nascent 
black hole. Since the collimation of the outflow in 
the double neutron star merger scenario should be 
mainl y contributed by the interaction with accretion 
torus (Al ov et al. H 2005"). a small Rq may be necessary 
to be consistent with the very small opening angle of the 
ejecta 9j found in the afterglow modeling. 

The GRB ejecta is mainly launched through the 
pole region of the rotating black hole. Clearly the out- 
flow should be from a site above the horizon surface. 
Along the pole that means i?o > ''g = '^GMbh/c^ ~ 
3 km (Mbh/1 Mq), regardless of the unknown spin of 
the black hole, where r^ is the Schwarzschild radius and 
G is the gravitational constant. In reality Rq is larger 
than Tg by a fact or of a few, as found in the numerical 
simulations fe.g iPopham et al.1 119991 : lOu et al. II2007I : 
IZalamea fc Beloborodov 1120111 ). With eq.^, for GRB 

090510 we then have Mbh < 22 M© Tpjf 3 gFi^^/i^, i.e., 
it is likely a stellar black hole at the center, in agreement 
with the compact-object merger model. 

For other short bursts, due to the lack of a robust 
estimate of Tint and then Fph, a reliable constraint on Rq 
is not possible. 

4. DISCUSSION AND CONCLUSION 

The central engine of Gamma-ray Bursts is widely 
believed to be one nascent stellar black hole surrounded 
by a hyper-accreting disk. With the current observa- 
tional data it is however not easy to pin down the physi- 
cal parameters, for example, the mass and the spin of the 
central black hole, the accretion rate (or alternatively the 
accretion-disk mass), and the initial radius of the outflow 
where it is luanched (Rq). The main reason is that the 
central engine hides deeply behind the electromagnetic- 
radiation surface and it is hard to break the degenera- 
cies between the parameters with the very limited ob- 
servational constraints. For some long bursts, in par- 
ticular GRB 090902B, people get a reliable estimate of 
Rq with the identified thermal spectrum component in 
the prompt spectrum. However for short events the lack 
of a reliable identification of such a component renders 
a reasonable estimate difficult. In this work we discuss 
whether we can estimate the disk mass in the specific 
scenario of binary-neutron-star merger (see section [2] for 
details). We also outline how to constrain Rq without 
the identification of an ideal thermal spectrum compo- 
nent and then applies it to GRB 090510 (see section |3] 
for details). Our main conclusions are the following: 

• Our semi-analytical estimate suggests that the 
nascent black hole formed in binary-neutron-star 
merger scenario rotates very quickly and the spin 
parameter is insensitive to the initial mass ratio of 



of these data sugge sts a typical Rq ~ 10^ cm IIRvde et al. II2OIOI : 
IZhang et al. Il2011h . We speculate that the intrinsic outflow was 
launched at a much smaller radius and was "choked" by some ma- 
terial at > 10^ cm, (i.e., the fireball was re-bom, see also Thompson 
et al. 2007; Ghisellini et al. 2007; Lazzati et al. 2009) for the fol- 
lowing two reasons. One is that the cooling of the disk material 
is dominated by neutrino radiation process, crucial for launching a 
baryonic outflow, only inside a radius ~ 10^ cm (jNaravan et al. I 
120011) . The other is that at a radius ~ 10^ cm, the annihilation 
of the neutrinos and anti-neutrinos should be very inefficient since 
the number density of (anti-)neutrinos drops with radius sharply. 



the double neutron stars (see section [52]), in agree- 
ment with t he results of the rec ent numerical simu- 
lation (e.g., iKtuchT^F^OIlOO^. Together with the 
finding that there is no significant mass ejecti on in 
the merger process (e.g.. lRosswog et al. 1 119991. the 
mass of the formed black hole M^^ as well as the 
spin parameter a may be reasonably deduced. Con- 
sequently a rough estimate of the accretion-disk 
mass is possible. 



As found in section 12.31 for about half of short 
GRBs in our sample, the disk mass is estimated 
to be ~ 0.01 — 0.1 Mq, in agreement with that 
found in the numerical simulation of the merger of 
binary-ncutron-star. For some other bursts, such 
as GRB 051221A and GRB 050724, a massive disk 
(^ several 0.1 Mq) is needed. This puzzle can be 
solved if the outflows of these short GRBs were 
launched via the more efficient magnetic processes 
rather than the neutrino mechanism or alterna- 
tively these short GRBs were from the neutron 
star-black hole merger (for which a massive disk is 
plausible). Since no significant divergence between 
the disk mass inferred from the observational data 
and that obtained in the numerical simulation has 
been found, we suggest that the compact object 
merger scenario for a good fraction of short bursts 
is viable and these events are promising gravita- 
tional wave radiation sources. 

• For GRB 090510, the initial radius of the outfiow is 
estimated to be < 6.5 x 10'^(rph/2000)^'* cm. Such 
a small Rq suggests that the neutrino- antineutrino 
annihilation products were mainly deposited in a 
small cavity surrounding the nascent black hole, 
as expected. Moreover, the small Rq imposes a 
constraint on the mass of the central engine M-^^ < 
22(Fph/2000)~'* Mq, consistent with the compact- 
object merger model. 

With the future short-burst-associated gravitational 
wave data, the binary-neutron-star merger model will be 
directly tested. Moreover the formation process of the 
disk, total mass and the mass ratio of the double neu- 
tron stars involved in the merger a nd the mass of the 
formed disk can be well constr ained (jKiuchi et al. II20T0I : 
iKobavashi fc Meszaros I [20031 ) . Consequently the valid- 
ity of our simple approach outlined in section [2] will be 
unambiguously tested. Considering that the magnetic 
process is usually much more efficient than the neutron- 
antineutrino annihilation to extract the energy, comb- 
ing the derived a, M^^ and Afdisk with -Eout, the nature 
of the outfiow-launching-process (magnetic or neutrino- 
antineutrino annihilation) will be reliably probed, too. 
For example, if the Mdisk inferred from the gravitational 
wave data is so small that can not produce the observed 
burst via neutrino mechanism, the magnetic outfiow- 
launching-process will be favored. 

Finally we'd like to point out that for the possi- 
ble short event G RB 100816A, t he pre liminary afterglow 
data reported in lOates et al. I (|2010[ ) tentatively favor 
the free- wind medium model, in which the progenitor 
should be a massive star rather than a pair of compact 
objects. Careful analysis of available optical/infrared af- 
terglow data is thus encouraged. If our speculation has 



been confirmed, the coUapsar origin of GRB 100816A 
with an intrinsic duration '^ 1.4 s would be firmly estab- 
lished, in support of the hypothesis that collapsar can 
produce short events and the progenitors of short GRBs 
are diverse. 
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TABLE 1 
Short GRBs and the mass of the accretion disks. 



GRB 


T9o{s) 


z 


■Et.Iso (erg) 


-Ek.iso (erg) 


6lj (rad) 


Afdisk (MQr 


References 


050509B 


0.05 


0.2248 


4.5 X lO"** 






0.02 


1 


050709 


0.07 


0.16 


6.9 X lO-is 


3.7 X 10^0 ° 


0.21 


0.03 


1-2 


050724 


3 


0.257 


4 X 10^0 


6.0 X 10^0 " 


0.2 


0.37 


1-2 


051221A 


1.40 


0.5465 


2.4 X 10^1 


1.0 X 1052 6 


0.10 


0.46 


3-4 


061006 


0.4 


0.4377 


2.1 X 10^1 


1.8 X 10^1 = 


~0.11 


0.10 


5-6 


070714B 


3 


0.9224 


1.2 X 10^1 


3.6 X 10^1 = 


> 0.08 "^ 


> 0.24 


7 


071227 


1.8 


0.381 


5.8 X 10^^° 


5 X 10^° " 


> 0.09 "^ 


> 0.13 


8 


090426 


1.25 


2.609 


3 X 10^1 


8.7 X 10^2 


0.06 <* 


0.30 


9 


090510 


0.30 


0.903 


1.2 X 10^^ 


5 X 10^53 


0.006 


0.06 


10-12 


100816A 


2.8 


0.8035 


5.8 X 10^1 


1.1 X 10^2 c 


> 0.01 '^ 


> 0.07 


13 



° We take the lower value obtained in lPanaitescu I II2006I ). 

' There was a flat segment in the X-ray afterglow and its origin is still unclear. Here we take an -Ek,iso ~ 10^^ erg required in the 

two-component jet model l lJin et al. 1120071 ). 

"^ This parameter is estimated by eq.l|4)l by taking Ec ~ 0.1, ea ~ 0.01 and Y ~ 0{1). 

'^ The jet opening angle is estimated by eq.ljSj- For GRB 090426, a jet break time t = 0.4 day (Nicuesa Guelbenzu et al. 2011) and 

n ~ 10 cm~^ (Xin et al. 2011) have been adopted. For other events, we take the time of the last Swift XRT detection, reported at 



to set a lower limit on the half opening angle. 
: 0.78, T -- 



|http://www.swift.ac.uk/xrt — curves/ as tj and n ~ 0.01 cm 

"^ We estimate Mdiakwith eq.JS}, adopting J\/bh = 2.7 Mq ^ 

— References: (1) Ex et al. 2005; (2) Panaitcscu 20 (3; (3) ISoderberg eFal . 
(&) IGol enctskii et al. ' 2006: (7) Ccn ko et al. Il200a f8l iCaito et al. Il2010l : (%) fXin et 
120111 ; (12) Gao ct al. 2009 ; (13) Gates et al. 112010 . 
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